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Abstract 

We derive a general formalism with which it is possible to obtain the time (r) 
dependence of the echo size for a spin in a stochastic field environment. Our 
model is based on "strong collisions" . We examine in detail three cases where: 
(I) the local field is iwo, (H) the field distribution is continuous and has a 
finite second moment, and (III) the distribution is Lorentzian. The first two 
cases show a T2 minimum effect and are exponential in as r ^ 0. The last 
case can be approximated by the phenomenological expression exp(— [2r/T2]^) 
with 1 < /3 < 2, where in the r — > limit (3 = 2. 
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Spin echo decay (SED) measurements, also known as T2, are conducted by a variety of 
experimental techniques, such as RF-/iSR 0, ESR 0, NQR, and NMR 0. With the recent 
explosion of high-Tc superconductivity research, NMR-T2 measurements in particular are 
receiving renewed attention, since they are very successful in probing both the normal 0] 
and superconducting states of cuprates. These experiments lead to a revival of theoretical 
activity, focusing on the calculation of the SED waveform for different sources of interac- 
tions such as spin lattice coupling, spin-spin coupling, and stochastic fluctuations. For this 
purpose, a variety of analytical and numerical models were applied. However, several 
dynamical features, observed experimentally, have not been accounted for. In this paper 
we provide new insight into these features by re-examining the echo decay waveform of a 
spin in a stochastic field environment, and use an analytical approach based on the "strong 
collision" model (see below) to yield quantitative understanding of SED. 

An earlier exact treatment of the stochastic problem, based on a diffusion like model, 
was presented by Klauder and Anderson (KA) . They found that for Lorentzian diffusion 
the waveform is Gaussian, and for Gaussian diffusion the waveform is exponential in (see 
bellow). Although the KA approach is physically more intuitive, the final result lacks three 
features: (I) the waveform does not depend on the diffusion rate, and therefore it cannot 
change continuously (for example, as a function of temperature), (II) they could not account 
for stretched exponential relaxation exp(— [2r/T2]'^) with (3 < 2, and (III) the SED depends 
monotonically on the diffusion rate, although it is natural to expect that when the diffusion 
is either very fast or extremely slow, the echo does not decay. As we shall see, our derivation 
allows for all these phenomena, and therefore might be applicable to some cases to which 
the diffusion model is not. In addition, we examine one case which was not considered by 
KA, namely, the local field is ±ujq. For this case our derivation is exact. 

In echo NMR, NQR, and ESR transverse relaxation measurements, a 7r/2 pulse is applied 
to a system of spins polarized along the z direction. As a result, a net polarization along 
the X direction (M^.) is obtained. In RF-/iSR the muons enter the sample with their spin 
already polarized along the x direction. After the pulse (or muon arrival), the spins evolve 

2 



with time, each one in its local field Bz, until time r when a vr pulse is applied, sending the 
X component of each spin 5*^. to —5*^ (and to —Sz). The spins then continue to evolve, 
and if Bz is static, an echo is formed at time 2r. If, however, the local field is dynamic, the 
phase acquired by the spin before the vr pulse is not necessarily equal to the phase lost after 
it, and the echo size diminishes as a function of r. This situation can be quantified by 



Ma;(2r) = M^.(0) ^cos j\{t)dt - j^^ uo{t)dt 



2t 

(1) 



where uit) = ■yBz(t), 7 is the spin's gyromagnetic ratio, and () is an average over all possible 
frequency trajectories. 

First we would like to evaluate Eq. |l| to lowest order in r. Assuming that the argument 
of the cosine is small, we can expand it to second order, and then evaluate terms such as 
Jq £ dt'dt" {uj{t')uj{t")) and /J" /^^^ dt'dt" {uj{t')u{t")) . Assuming a correlation function of the 
form 

{uo{t')uo{t")) = (uj^) exp (-Z/ \t" -t'\) = (uj^) \t" -t'\ + ...), (2) 

where (u;^) is the second moment of the instantaneous field distribution, we find 

M,(2r) = M,(0) (^1 - ^ (iu^) uT' + ..}j. (3) 

Equation ^ is well known and will serve as a test of our derivation. 

Next we shall evaluate Eq. |l| to all orders in r by making some assumptions concerning 
uj{t). We quantify the dynamical fiuctuation using "indirect echo" and the strong collision 
model. Indirect echo is equivalent to the situation described by Eq. |l| but instead of Sx —>■ 
—Sx at the TT pulse, the frequency is reversed {00 —uj); in Fig. |l|a we demonstrate indirect 
echo by showing that a reversal of u; at r leads to Sx{2t) = ^^^(O) = 1. The strong collision 
model accounts for uj(t) by allowing frequency changes only at specific times ti,t2 . . .tn- The 
probability density of finding the frequency uj at any time interval is taken to be the line 
shape p{uj). A demonstration of this situation for a particular spin is presented in Fig. |^b. 
Here the spin has experienced two frequency changes at times ti and ^2 before the vr pulse 



and one change after the n pulse at t^. As a result Sx{2t) ^ SxiS^^ and on average the echo 
size will decrease as a function of r. This type of dynamical process results in a correlation 
function in the form of Eq. 0. By comparison, in KA's model the frequency after each change 
depends on the frequency before the change. 

We shall now treat the case of an ensemble of spins and average over all possible field 
changes, the times at which they take place, and all possible fields in each time interval. If 
there are n hops at times ti^ . . . ,tn before the vr pulse and m hops at times t„+i, . . . , tm+n 
between the tc and the observation time t = 2r, the phase acquired by the spin {9n,m) is 

m 

Gn,m = ^n+m+l{t ^ tn+m) + '^^^j+niin+j ~ tn+j-l) (4) 

i=2 



- cc;„+i(t„+i - r) + co'„+i(r - t„) + ^uji{ti - 

1=1 

The polarization along the x axis is therefore 

where 3fJ stands for the real part; we shall omit it from now on. We first average over all 
possible frequencies uji in the time segment [ti-i,ti] and define 



Mx{t, T]ti, . . . , tn+m n+m+1 

)duJn+m+lMx{uJi, . . . 5 ^^n+m+l] t,T]ti, . . . , tn+m+l)- 

This results in 



Mx{t, T]ti, . . . , tn+m) — git — tn+m) 



n gitj+n — tj+n~l] 
i=2 



g{2T — tn+l—tn-l] 



\{g{ti -ti_i) 



.i=l 



where g{t), also known as the free induction decay (FID) function, is given by 

roo 

g{t) = / p{uj) exp{iujt)dLU. 



The probability density of finding exactly n + m hops at times ti, . . . , tn+m is 

n+m n+m 



exp [-z/(t - tn+m)] n [-^i^i " ^^ti = z/"+'" exp(-z/t) Yl dti, 

i=l i=l 

where u is the field selection rate. Thus, the averaged spin polarization at time t is given by 



oo oo 



M.it) = E E ^^"+™exp(-z/t)/„,„(t,r), (5) 



71=0 m=0 

where 



/ dtn--- dtiM^{t,T;ti, . . . ,tn+m)- (6) 

Jt Jt jo jo 

The integration hmits guarantee that tj+i > t,. 

We can simphfy Eq. |^ by turning the time at which the vr pulse is apphed (r) into a 
running variable (t') whose value is fixed with a 6 function. The 6 function should force the 
sum of time segments from zero until r to be equal to the sum of time segments from the r 
until 2r, namely, 

(n m+1 
(f - 1„) + E(ii - u-i) - E (Wi - tn+j-i) - (Wi - 1') 
i=l j=2 

where tn+m+i stands for 2r. As a result 

r^T rtn+2 rtn+i rt' rt2 

In,m{'^T,T) = C?t„+„ ■ ■ ■ / dt^+l dt'l dtn,--- dtiM^{2T, t' ] ti . . . tn+m)S{t' - t) 

JO JO JO JO JO 

(7) 

and the integrand in Eq. |^ is a function of time differences only. 
We now introduce the integral representation of the 6 function 

1 

S(x) = — / exp{inx)dn, (8) 
and the Laplace transform of M^: 



oo 



MJs) = 2 M^(2r)exp(-2sr)rfr. (9) 

JO 

By inserting Eq. ^ into Eq. ^ Eq. ^ into Eq. ^ and substituting this in Eq. |^ we find that 
all the integrals decouple and 

1 /"OO ^ ^ 

M.{s) = - dnUz.,z^)Y.T.i^hi^-)Ti^hi^^)T m 

Z,Tl J —oo „ n ™ n 



n=0 m=0 



where 



z± = s + u ± iQ/2, 



(11) 



poo 

fi{z±)= duexp{-z±u) g{u), (12) 
Jo 

and 

;,(,_,,j = W£z)±il(£±). (13) 

Finally, |z//(z)| < 1, and performing the sum in Eq. ^ gives 

=^r — .f'^tuf^^ ,1 (14) 

from which we obtain the time dependent nuclear magnetization by 

M,(2r)=/:-i(M,.(s))^^^^, (15) 



where is the inverse Laplace transform operator. Using Eqs. 11 to 15 one can obtain 
the echo decay knowing only the line shape or the FID function. Now let us examine three 
simple cases: 

Ising field - this case materializes when the observed spin is coupled only to one unob- 
served spin 1/2 which fluctuates stochastically, leading to either field up or down at the 
observed site. In this case the field distribution is given by 

1 1 

and its second moment is (cu^) = uJq. Here we have to be careful about the definition of u 
since the field selection rate z/, which appears in Eq. 0, is twice the rate u-j- at which there 
are actual field changes. It is that counts in our derivation, and we find 

Ml{2r) = (o;o^ + u^fi sin(2/ir) - ^ cos(2/ir)) , (16) 

where ff = uJq — u^, and the superscript I stands for Ising. For < the result is the same 
except that fi ^ i\fi\- The expansion of Eq. |TB| to lowest order in r agrees with Eq. ^ In 
Fig. Ha we present on a semi-log scale as a function of 2ujqt for various values of zz/wq. 
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It is clear from this figure that when either v/ojq <^ 1 or v/ojq ^ 1 the echo decay rate is 
weak compared to zz/cuq — 1- To quantify this phenomena we define T2 as the time at which 
the echo size decreases to 1/e. In the Ising case we find that T2 reaches its minimal value of 
3.146/^^0 when v/ujq = 0.69 . 

In the inset of Fig. ^ we depict M^. vs. (2u;ot)^ for v/loq = 0.69. This presentation 
emphasizes a surprising fact that when u/ujq ^ 1 the waveform looks Gaussian over 2 orders 
of magnitude in echo size, even though at early time it is exponential in r^. When fitting 
this Gaussian to M^(2r) = exp{-{2uJoT/T2G)^ /2) we find T2G = 2.17. 

A distribution with a second moment - It is useful to examine a continuous distribution 
with a finite second moment so as to compare with Eq. ^. One such distribution is: 

and its second moment is given by (cj^) = This leads to 

^'-(^r) ^ - ^^^^ (IS) 



4(ct-2,.)/J > ' 2f,(a-2v) 

where = (cr^ — 2(tz/ — z^^)/4. Again, for < the result is the same, except that 
fa ^ i\fa\- An expansion of Eq. |T8| around r = agrees with Eq. |^, thus demonstrating the 
validity of our derivation once again. In Fig. 0d we depict Eq. |l8]for various values of u/a. 



It is clear that the echo decay is strongest for u = 0.88cr ~ y (u;^) . In fact, at this value of 
z/, T2 is minimal and equals 5.75/cr. 

Lorentzian distribution - in this case the equilibrium distribution is taken to be 

and we find 

, , Aexp(— 2z/r) — z/exp(— 2Ar) , ^ 

mI;{2t) = — ^ (20) 

where L stands for Lorentzian. This expression has interesting properties. An expansion of 
Eq. ^ around r = gives 



M^^{2T) = l-\\u{2Tf + 0{T') 

which means that at early enough times the relaxation shape is Gaussian. One should note 
that this expansion does not contradict Eq. ^ since a Lorentzian does not have a second 
moment. However, for A ^ z/ the relaxation is exponential for Ar ^ 1 with the relaxation 
rate v. Similarly, when \ <^ v the relaxation is exponential for z/r ^ 1 with the relaxation 
rate A. This suggests that experimental data which stem from Eq. ^ can be well fitted to 
a stretched exponential 

r2rl^\ 



M^(2r) = exp 



(21) 



with 1 < /9 < 2. In Fig. ^ we depict three data sets of M^{2t) obtained from Eq. ^ for 
various values of u/X. Unlike in the previous cases, the Lorentzian case shows a continuous 
increase in relaxation rate with increasing u. In this figure we also demonstrate the best fit 
of the data sets to Eq. The fits are quite good over more than an order of magnitude in 
echo size, and when experimental data are fitted, Eq. ^ can easily be confused with Eq. |21|. 
In Fig. ^ we show the parameters /3, and 1/(AT2) as a function of z//A. While T2 decreases 
monotonically with increasing fluctuation rate, the power (3 goes through a maximum at 
Z//A = 1. However, it should be mentioned that the value of (3 depends on the range which 
is used for the fit. 

It is interesting to compare our Lorentzian result with that of KA. In the KA model 
the field dynamics at the site of the observed nuclei is generated by flipping some other 
unobserved individual spins. Therefore, in their model, it is more likely to undergo small 
field changes than large ones. The situation KA tried to describe could still be approximated 
by the strong collision model if z/ ^ A since then many unobserved spins are flipped before 
the observed nuclei evolve considerably with time. This suggests that in reality, for z/ ^ A, 
we should expect /5 = 1, as found here, and for z/ ^ A we should expect /3 = 2 as found by 
KA. Between these two limits (3 should change continuously. 

We thus provide a recipe for obtaining the time dependence of the echo size for a given 
field distribution. We examined three particular cases and found a natural explanation 
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for experimental and conceptual features, such as stretched exponential relaxation and T2 
minima, which have not been explained quantitatively before. 
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FIGURES 

FIG. 1. A demonstration of the indirect echo when there are no local field fluctuations (a), and 
when the field is dynamical (b) and changes instantaneously. 



FIG. 2. The echo decay for: (a) the Ising field (Eq. pig ) plotted against 2ujqt, (b) A continuouse 
distribution with a second moment (Eq. 18) vs. 2aT, and (c) Lorentzian field distribution (Eq. |2C 



as a function of 2Ar. The solid line in panel (c) represents a fit to Eq. 21 as described in the text. 



FIG. 3. The parameters 1/(AT2) and (3 which allow the best approximation of Eq. 2C with 



Eq. 21. In the inset, the Echo decay for the Lorentzian case with u/X = 0.7 is shown. 
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